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ABSTRACT 


The principal component analysis enhances 
the contrast existing between the different cover 
types present In an Imagery. A procedure Is prese- 
nted In this paper with regards to the determina- 
tion of the principal components. The method Is 
tested for a portion of the LANDSAT Imagery perta- 
ining to Anantapur region. Another technique, using 
the concept of non-linear contrast stretching is 
defined and developed and carried out on the same 
Imagery. The results are presented as photographs. 
An Interpretation of the geology of the region Is 
derived from these photographs. 


1. IHTHDDUCnol 

Tbr processing of LANDSAT Imagery usually three distinct techniques are 
used. One technique Is to process the Imagery with the help of computers, the 
second one is to Interpret the imagery visually, and the third is a hybrid 
approach where both autematlc and visual Interpretation are used. 

In computer processing of LANDSAT MSS Imagery, there are two steps. The 
first step is to preprocess the image in order to remove noise, and to Improve 
contrast. The second step Is to use the multi spectral Imagery and classify the 
cover types. The multlspectral classification is very effective for some speci- 
fic applications like study of landuse, forestry and agricultural inventory. 

By this classification It Is possible to generate thematic maps also (1). 

In some applications like In geology, thematic mapping Is not of relevance 
because geological interpretation does not wholly depend on surflclal proper- 
ties and signatures as there Is a limited correlation between vegetative and 
other surface cover and geology. Fbr this reason the Interpreter or the photo- 
geologist usually prefers to study the Image Itself visually in order not to 
leave out any significant data or clue. However In some cases vlsvial Interpre- 
tation of the Imagery Is not possible because of the low contrast existing bet- 
ween different objects present In the Imagery or due to some Inherent drawbacks 
of the photographic recording devices. To avoid these drawbacks various numeri- 
cal techniques have been used to enhance the satellite acquired data. Some of 
the commonly used tediniques are contrast stretching, band-to-band ratlolng ^md 
spatial filtering. 

In this paper an attempt Is made to develop procedures for enhancing the 
contrast between different objects In MSS Imagery. These procedures vlU be 
useful for the analysis and Interpretation of L^TDSAT Imagery for geological 
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applications. As a result of these enhancement techniques, one can obtain cert- 
ain information which is otherwise not available for visual interpretation. 

The basic idea of one of these procedures is the use of the concept of 
principal component euialysis as discussed in (2,3). The method allows for shift- 
ing of gray values in the photograph after getting principal components. The 
analysis depends on the multispectral classification approach but additionally 
processes the signals so as to achieve better contrast in the picture. 

The second procedure described in this paper also uses the multispectral 
classification technique and evolves a non-linear weighting function to obtain 
contrast stretching. 


2. PROCKDUBE FOR AHALYSIS 

The method of analysis can be considered in three steps namely the classi- 
fication of the nniltispectral data, the computation of covariance matrices and 
finally the determination of principal components. 

2.1. CLASSIFICATION 

The first step is to classify the area considered, using the available 
multispectral data. The classification procedure will identify and separate out 
the different cover types. An unsupervised classification technique is used for 
this purpose (4). 

2.2. COVARIANCE MATRICES 

After completing the classification, the next step will be the computation 
of the covariance matrices for the different cover types. Because of inherent 
spectral correlation which the surfaeial features of theaarth possess in diffe- 
rent spectral channels, the covarieinee matrices will have invariably non-zero 
elements. The covariance matrix is denoted by CL where 'm* represents cover 
type. 


2.3. PRINCIPAL COMPONENTS 

The last step will be the determination of Principal Components. This step 
involves the oomputation of eigen values and the corresponding eigen vectors Ibr 
each cover type. The eigen values are arranged in the descending order. 

The principal components are obtained from the original N spectral bands 
by the following transformation 

Y \ ( X-U„) (1) 

where X is the N-dimenslonal vector of spectrml intensities associated 
with each pixel or ground resolution element, 

Uju is the mean vector of dimension N for a cover type *m', 

Y is the vector of principal components, and 

TL is the N X N orthogonal matrix derived from the covariance 
matrix C^ of each cover type 'm*. 

The rows of ^ are the N normalised eigen vectors of . 

Using Bqn. (1), new gray values for all the pixels are obtained. The pro- 
ced\ire will generate as many number of principal components as there are number 
of spectral channels. Essentially a new set of pictures are generated from the 
given four cheuinels by mapping Into new four components. 

The covariance matrix of the prlnclp€il comjxjnents is then Gy = T^ CL T^ = 

[i-PJ 
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where are the rarlances of the principal component and are also 

the eigen values of C • It is evident from C that the prlndpeil components 
are vincorrelated and “ each component has variance less than the previous 
component* 

2.4. NOH-LINEAR CONTRAST STRETCHING 

In order to enhance the contrast of Imagery besides the principal compone- 
nt technique, various techniques have been tried. Bond and Thomas (5) have 
reported a method that performs contrast stretching. In a similar fashion, but 
using a different criterion, a technique for contrast stretching is formulated 
here. In this new gray values for pixels belonging to each cover type are obta- 
ined using non-linear stretching over the gray values of the original channels 
of LANDSAT imagery. 

Using the knowledge that a particular pixel belongs to a certain cover 
type, the gray value of this pixel is modified to a new value, by a weighting 
factor. This weighting factor, is given by 


G (j, k) = 


P (J,k) 

t P(J.k) 

k=1 


(3) 


where G ( j ,k) 
P (J,k) 


is the weighting factor for 'k*th class in the 'J’th band 
or channel, and 

is the mean value of reflectance of the’k'th cover type 
in the ' j 'th band. 


3. IMPLEMENTATION OF THE PROCEDURE 

The analysis procedures as described in the previous section are carried 
out on LANDSAT imagery, pertaining to Anantapur region (I. D. No. 2053-04265, 
l6th March 1975). 

3.1. SITE FOR STUDY 

A portion of the imagery covering an area of 30 Rn x 30 IQn is considered 
for this purpose. The area includes rocks of varied lithological entities like 
granites, qmrtzltes, limestones 2 uid shales, shales and sandstones, and basic 
igneous intruslves. Penner is the major river flowing through the test site. 
Central belt of the area trending NNW-SSE comprises hills of varying heights 
and different rock units. Eastern part of the area is almost plain and shows 
the presence of basic dykes. 

Figures 1 and 2 represent gray values of LANDSAT bands MSS 5 and 7 respe- 
ctively. Sinuous course of the river Penner can easily be deciphered in Figure 
2 than in Figure 1 due to its contrast in gray tone from the surroundings. 
Central belt of hilly terrain corresponds to limestones and shales and appears 
more clear in Figure 2. The hill range at the north-western part of the area 
is perhaps less ^slble in gray value photographs of bands 5 and ?• 

3.2. DATA ANALYSIS 

"Bie portion of the imagery has been digitized using the Optronics P-1000 
photoscan system. The digitizer uses 256 gray levels, gives out one byte of 
data for each resolution element. !Die digitized data for each spectral band is 
stored in a sequence on a CCT. For convenience of handling 2 uid analysis, the 
data is transferred to a magnetic disc 

Fbr analysis purposes, the registration of the imagery in all the b 2 mds is 
an important step. The Sequential Similarity Detection Algorithm (4,6) is used 
to achieve the correct registration. The analysis of the imagery is performed 
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In segments, each segment covering an area of 10 Km x 10 Km having ^0,000 
pixels. Each pixel will have four components corresponding to the s ctral 
intensities in the four channels. 

3.2.1. THE PRINCIPAL COMPONENTS. Bie inherent differences existing bet- 
ween different cover types is utilized in the derivation of principal componen- 
ts. 


The data present In the above area is classified using the unsupervised 
classification algorithm (6). Tte analysis is carried out in nine segments each 
of 100 sqxmire kilometer. In a typioil segment there are six cover types, the 
statistics of which is presented in Thble I. Die covariance matrices of eill the 
cover types are then computed. The computed covariance matrix for a cover type 
is given in 'Gable II. 


TABLE I. STATISTICS OF A TYPICAL SEGMENT 


Cover 

Type 

Popula- 

tion. 

Mean 

Reflectance 

St2usdard 

Deviation 

kSS 4- 

MSS ^ 


MSS 7 

MSS 4- MSS 5 MSS 6 



1827 

220.7 

192.2 

183.6 

206.1 

7.7 

10.2 

10.9 

8.7 


2981 

214-.7 

173.8 

165.0 

194.9 

7.5 

7.6 

9.0 

7.5 


2401 

207«4 

15^.7 

152.0 

185.8 

8.4 

7.4 

8.2 

6.9 

■■ 

1405 

203.9 

152.0 

139.5 

171-3 

9.1 

10.3 

8.4 

7.2 


1098 

217.6 

187.6 

153.2 

178.0 

7.9 

13.6 

8.6 

9.2 

H 

288 

217.4- 

196.4 

170.4 

187.9 

20.7 

37.7 

31.1 

26.3 


TABLE II. COVARIANCE MATRIX, EIGEN VALDES, AND 
EIGEN VECTOR MATRIX FOR A COVER TYPE. 


COVARIANCE MATRIX 

EIGEN 

VALDES 

EIGEN VECTOR MATRIX 


'56.48 7.67 -7.61 -2.GO" 

7.67 58.72 -3.37 1.52 

-7.61 -3.37 80.54 -3.74 

-2.00 1.52 -3.74 56.28 


84,32 

62.64 

56.82 

48.23 


'-O.31 0.56 -0.21 0.74' 

-0.22 0.73 0.35 -0.54 

0.92 0.33 0.12 0.16 

_-0.11 -0.20 0.90 0.36 



The new gray \^lues for all the pixels are then detennlned using the 
matrix t ran s forme tl on given by Eqn. (1). !Sie principal components are then 
obtained 2Uid are stored on a magnetic disc in a format of a total segment. 

These are later used for writing the image on a film transparency in each of 
the principal components. 

Figures 3 and h show the prints of first two principal components obtained. 

3.2.2. CONTRAST STRETCHING. As explained earlier, an alternate scheme 
for contrast stretching is to use the information regarding the cover tyi>es to 
%dileb the pixels belong. Ebr this purpose, the gray values of the contrast 
stretched imagery are obtained by multiplying the original gray values of the 
pixels by the weighting factor, defined by Eqn. (3)* 
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Figures 5 and 6 show the prints of contrast stretched imagery 1 and 2 
respectively. 


h, VISUAL INTERPRETATION FOR GEOLOGICAL STUDIES 

The visual Interpretation discussed here is supported by limited field 
checks of the area under study. Three major units, namely hills comprising 
sedimentary rocks, water bodies and daiic gray belts corresponding to dykes, are 
considered for assessing the Interpretabllity of MSS bands 5 and 7, principal 
components 1 and 2 and contrast stretched Imagery 1 and 2. 

‘f.l. INTERPRETATION OF PRINCIPAL COMPONMTS 

Using the photographic data from Figures 3 and V corresponding to principal 
components 1 and 2 respectively, a visual interpretation of the geology of the 
area is carried out. Central belt of limestones and shales is represented in 
almost vmiform gray tone in principal components 1 and 2. River Penner is dis- 
cernible in principal component 1. Due to the sharp tonal contrast the hill 
range at the north-western part can be picked up easily in principal components 
1 and 2 in comparison to gray value pictures of MSS 5 and 7. 

Figure 7 shows the interpretation. A band of quartzite separating limesto- 
nes and shales from granites is delineated in Figure 7« In the eastern part of 
the area, dykes running NW-SE are detected in shales. Same can also be seen in 
principal component 2 where it is less visible. Due to very little gray tone 
variations, it has not been possible to delineate these units from principal 
components 3 and 4. For this reason the prints of the nrindpal component 3 and 
4 are not Included in this paper. 


4.2 INTERPRETATION OF CONTRAST STRETCHED IMAGERY 

River Penner and the hill range at NW comer appears more clearly In 
Figure 6. In addition, several hill and sedimentary rocks can be easily picked 
up on Figures 5 and 6 of the contrast stretched Imagery. 

Using the photographic data from Figures 5 and 6, a visual interpretation 
map is generated and is shown in Figure o. A quartzite hill range separating 
limestones and shales from sandstones and shales, and basic Igneous intruslves 
are clearly discernible on the enhanced images. These are delineated in Figure 

8 * 


5. CONCLUSION 

Two methods are presented to enhance the contrast among different cover 
types present In an Imagery. It is noted that the principal component analysis 
has compacted all the information content present in the Imagery especially to 
the first two components. The non-linear contrast stretching method has signi- 
ficantly improved the contrast among the cover types of the area under study. 
Certain geological formation like quartzite, limestones 2 uid shales, and basic 
igneous intruslves are more clearly seen in the principal components and also 
in the contrast stretched Imagery. 
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FIGDRE 1. GRAY VALDES OF MSS 5. 


FIGURE 2. GRAY VALDES OF MSS ?. 
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FIGURE 5 * CONTRAST STRETCHED IMAGERY 1 


FIGURE 6 


CONTRAST STRETCHED IMAGERY 2 
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Fijur*?- VISUAL INTERPRETATION FROM PRINCIPAL COMPONENTS 



Fijure 8 -VISUAL INTERPRETATION FROM CONTRAST STRETCHED 

imagery 




